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Abstract 

We show a strong indication of the existence of a large Fermi surface in 
the one-dimensional Kondo lattice model. The characteristic wave vector 
of the model is found to be kp = {1 + p)vr/2, p being the density of the 
conduction electrons. This result is at first obtained for a variant of the 
model that includes an antiferromagnetic Heisenberg interaction Jh between 
the local moments. It is then directly observed in the conventional Kondo 
lattice [Jh = 0), in the narrow range of Kondo couplings where the long 
distance properties of the model are numerically accessible. 
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I. INTRODUCTION 



Many experiments tell us that one of the low temperature states of the heavy fermion 
materials is a Fermi liquid whose quasiparticle masses are 10^ to 10'^ larger than those of the 
normal metals [Q. It is commonly believed that the heavy Fermi liquid state is one of the 
possible ground states of the Kondo Lattice model (KLM). The most popular description 
of the heavy fermion state is that below a characteristic temperature Tcoh (Tcoh ^ lO-ft'), 
the conduction electrons (c— electrons) and the local moments (/—electrons) have common 
excitations. This picture leads to the important fact that the localized spins also participate 
in the Fermi surface (FS). Thus, the FS has a large area. De Haas- Van Alphen measure- 
ments 1^ of some heavy fermion compounds have shown that the FS has a / character. 
Ab-initio Local Density approximation (LDA) computations of the band structure of 
heavy fermions also predict a large FS, although these calculations fail to reproduce heavy 
masses. These results are indeed expected from the Luttinger theorem which states that the 
volume of the FS is unchanged by the electron-electron interaction. Intuitively, however, it 
is not straigthforward to understand in the framework of the KLM, how the /-electrons can 
be included in the FS since there is no explicit hybridization between the c-electrons and 
them. LDA computations have shown the existence of a narrow / bandwith which is due 
to hybridization The KLM itself is an effective model of the periodic Anderson lattice 
(PAM) in the limit of nearly integral valence and of strong Coulomb correlation. One may 
thus wonder wether or not the residual /-electrons itinerant character present in the strong 
coupling limit of the PAM plays a crucial role in the formation of the heavy Fermi liquid 
state. 

Analytical [§] and numerical 01 PI studies have focused on the FS of the ID 
KLM. Their results remain very controversial. In the following, we will study the ID 
KLM numerically. A numerical study of the ID KLM presents two essential difficulties. 
The first one arises from the very low energy scale of the Kondo physics which requires the 
investigation of lattices of very large sizes. In real materials, the heavy masses involve a very 
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small value of the quasiparticle weight. For the ID model, in the physical range of parameters 
{Jk ^ 1); the size of the expected singularity in the electron momentum distribution nc{k) 
is likely to be very small. The second problem is the occurrence of a ground-state phase 
transition from a paramagnetic (PM) state at weak couplings to a ferromagnetic (FM) 
state at strong couplings ^ 0. Therefore, the results of the strong coupling regime where 
the model converges more rapidly to the thermodynamic limit cannot be extrapolated to 
the weak-coupling region where size effects are still significant even in very long chains. We 
will show that the study of a KLM (Eq. 1) in which the strong coupling regime is smoothly 
connected to the weak coupling one can give insight of the existence of a large FS for the 
usual KLM . Then, a careful analysis of the usual KLM will be made. A large FS means 
that the Fermi wave vector is located at kp = kp^ + ^, kp^ being the Fermi wave vector of 
the c— electrons only. We will consider the following KLM: 

H= -t^icf^Ci^^^ + h.c.) + JK^Si^ ■ + Jn^Sif ■ Si^^f (1) 

is it 

Where = \ Es,s' ^ts^'ss'^is' ^^"^ ^ic is a localized spin. The hopping integral t is set 
to 1. A direct Heisenberg exchange term between local moments Jh is introduced here. 
We have recently argued |^ that the occurence of a ferromagnetic phase transition in the 
KLM in the strong coupling region is due to the fact that the RKKY interaction becomes 
ineffective. One would thus expect that a sufficiently strong antiferromagnetic Heisenberg 
coupling between the local moments can stabilize a PM ground state. In the conventional 
KLM, this term is usually omitted because typical lattice parameters in the heavy fermion 
compound are 3.5 — 4A while the ionic radii of the / ions are less than lA, so that the 
overlap between the / orbitals is negligible. It should be noted that the strong Jk regime 
of the Hamiltonian (1) is relevant in the study of High-Tc materials. In this case, Jh is the 
superexchange interaction between copper ions. The double occupancy in the c— electron 
band is naturally suppressed by Jk, so that one need not include an explicit repulsion term 
between these electrons. We have investigated Hamiltonian (1) using the density matrix 
renormalization group (DMRG) method |jlO|. We have chosen an algorithm with open 



boundary conditions. We keep between 64 and 150 states in the two external blocks. These 
states are labelled by the 2;— component of the total spin S^. The ground state corresponds 
to the lowest state with = 0. The maximum truncation error is in the order of 10~^. 
Although we have reached = 60 sites, the longest distance in the calculation of the 
correlation function is L = 22. Because we have first built lattices of 20 sites before we start 
to calculate the correlation functions. This way, we minimize the end effects and density 
fluctuations that are larger in the early steps of the algorithm. 



II. RESULTS FOR THE KONDO-HEISENBERG CASE {Jh 7^ 0) 

The Hamiltonian (1) has a PM ground state with a Luttinger Fermi surface when Jh = 
0.5 in the strong regime. We have chosen our value for the Heisenberg coupling on the 
grounds that it is necessary that Jh exceeds the effective FM coupling J^ff between the / 
electrons to stabilize the PM ground state for all Jk- Our first choice was Jh = 0.1. For this 
value, we found that the ground state is PM in the weak Jk region, FM for intermediate 
couplings, then PM for strong Jk- We can estimate J^jf by using the results of the strong 



coupling expansion [0. JeffJR/t^ is approximatively 0.05, 0.2 and 0.25 respectively for 
the partial band-fillings p = 0.25,0.5 and 0.75. Then, knowing that the strong coupling 
region starts for Jk of order 1, one obtains the upper bound, J^/f^ — 0.25 for t = 1. Our 
value of 0.5 thus include a security factor. In Fig.l, we show nc{k) (the Fourier transform 
of (^cf^Cj^')) for the band-fillings p = 0.25,0.5 and 0.75 for Jk = 10. nc{k) at 0.25 and 
0.75 are affected by density fluctuations, since the band- filling is not constant during the 
DMRG iterations. Nevertheless, clean singularities are observed at kp = 0.6257r, 0.757r 
and 0.8757r. The magnetic structure factor of the localized electrons Sf{k) (the Fourier 
transform of <^5'jyS'|j)), shown in Fig.2, presents a maximum at 2kF = 0.757r, 0.57r and 0.257r 
respectively. These values correspond to 27r — 2kF, since 2kF is greater than vr. Clearly, 
neither the bare c— electron nor the bare /—electron signatures are detected. There are 
instead unique compound-particles propagating with a characteristic wave vector at kp. An 
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analogy can be made with the Hubbard model |T^. At half-filling, in the Jk = oo limit, all 
the conduction electrons form on-site singlets with the localized spins, so that the overall 
system is in a singlet state. The non half-filled cases correspond to the introduction of holes 
in the system. These holes which can hop from site to site are associated with the N — Nc 
unpaired /—electrons, Nc being the number of c— electrons. Obviously, double occupancy of 
holes is forbidden : one has a. U = oo Hubbard model of ph = 1 — p hole density. Depletion 



effects ||T3| are also observed in Fig.2: the reduction of the conduction electron density 
increases the tendency to magnetism. The maximum of Sf{k) increases when p decreases. 

Now we wish to discuss how the system evolves when Jk is reduced. We show in Fig. 3 
nc{k) at p = 0.5 for Jk = 10,8,6,4,3,2.5 and 2. The height of the singularity at kp 
decreases as Jk is reduced. Concurrently the drop at kp^ which, was negligibly small in the 
strong Jk case, increases. The c— electron character is progressively enhanced. The local 
spin-spin correlation which is = —.3748 ~ — |p at Jk = 100, is equal to —.366 at 

Jk = 10 and —.183 at Jk = 2. The deviation of this quantity from the perfect on-site singlet 
value — |p means that the singlet clouds have a spatial extension at lower Jk- At Jk = 2, 
short-range effects coexist along with the long-range behavior of the system. Clearly, it 
becomes hard to define the exact position of the FS. This can be better illustrated in 
Sf{k) (Fig. 4) at p = 0.5. A new peak appears at A; = vr at small Jk, signaling short-range 
antiferromagnetic correlations. At Jk = 1-5 (not shown here), the peak at 2kF is not seen. 
We can no longer observe the long-distance behavior of the model because of the finite-size 
effects (long correlation length; see the discussion below). Since there is no phase transition 
in the system, the weak and the strong Jk regimes are continuously connected. We thus 
believe that the FS is large even at smaller Jk- The above discussion is similar to the 
one made by Kotliar in the framework of the PAM |T^. The local singlets of the strong- 
coupling limit are obtained as the Kondo resonances are pulled out of the c— electron band 
by increasing Jk- The FS is conserved during this process. 
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III. RESULTS FOR THE CONVENTIONAL KLM {Jh = 0) 



We now discuss if the above conclusions can be extended to the PM phase of the con- 
ventional KLM. At first sight, one can argue that the conventional KLM is adiabatically 
reached by taking the limit Jh 0. Thus the conventional KLM may have a large FS in 
its PM phase. It should be noted from the above results and from the knowledge of the 
occurrence of a phase transition, that the range of Jk in which we can expect to detect a 
large FS with numerical methods in the conventional KLM is narrow. An estimation of 
this range can be obtained by examining the Kondo coherence length C,k of the one-inpurity 
problem, = vp/Tx, where vp is the Fermi velocity and Tk the Kondo temperature. The 
observation of the long range properties of the model is only possible at distances r ^ ^k- 
For r < finite size effects dominate, only short-range effects governed by the RKKY 
interaction will be observed. S0rensen and Affleck have recently made an accurate computa- 
tion of illl. Some typical values are = 1, 2, 4.85, 8 and 23 for Jk = 2.5, 2, 1.5, 1.25 and 
1 respectively. Clearly, for Jk ^ 1.25 the long-range behavior of the model is not attainable 
since the longest distance in our study is L = 22. Moreover, the FM transition occurs at 
Jk — 1.5 for p = 0.5 and Jk — 2.75 for p = 0.75. At low band-fillings, depletion effects 
enhance the FM instability. The PM phase boundary is shifted towards weak couplings 
where ^k becomes very large. Hence, at quater filling where there is no density fluctuations, 
this range is very narrow. Thus we have chosen to study the KLM at p = 0.75 in the 
range 1.25 ^ Jk ^ 2.5. In Ref. we have found that nc{k) displays a sharp drop at kp^ 
and Sf{k) presents a maximum at 2/cp^ in the small Jk regime. As Jk was increased, these 
features vanished before the phase transition was reached. But we were unable to draw a 
firm conclusion about the location of the FS. A more careful analysis will now show that 
we had observed a short-range effect governed by the RKKY interaction. Furthermore, we 
will identify the true long-range properties of the model which are not easily observable. 
In Fig. 5, we display Sf{k) at p = 0.75 for Jk = 1.25,1.5,1.75,2, and 2.5. Starting from 
Jk = 1.25, we can only detect the RKKY maximum at 2kF^ = 0.757r. For Jk = 1-5, the 
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maximum of Sf{k) is still located at the RKKY wave vector. But one can also observe a 
local maximum at the position of the large FS at 2kp = 0.257r . At = 1.75 and 2, the 
height of the RKKY maximum decreases. At the same time, the height of the maximum 
at 2kF increases. This trend is unambigously confirmed at Jk = 2.5, where the maximum 
at 2k p is now the highest. The height of these maxima however, is still smaller than the 
one at 2kF^ for Jk = 0.5 in ref [9], and are thus harder to detect. Finally, nc{k), shown in 
Fig. 6, corroborates the existence of the large FS in the KLM. nc{k) drops monotonously 
when k < kp = 0.8757r, shows a small plateau just before k = kp, and then drops abruptly 
at k ~ kp. The width of the plateau shrinks and then becomes undetectable when Jk is 
decreased. At the same time, the drop at kp^ is enhanced indicating that the short range 
effects are becoming dominant. This is consistent with the results for Sf{k). 



IV. CONCLUSION 

In summary, we have shown that the KLM with a direct exchange Heisenberg coupling 
has a large FS in the strong Kondo coupling limit. We have shown that this phase is 
continously connected to the weak Kondo coupling regime. As a consequence, the latter 
model has a large FS in the small coupling regime. The conventional KLM can continously 
be reached by taking the limit of vanishing Heisenberg coupling. We have concluded from 
this that its FS should have a large area in its PM phase. Direct numerical computations 
made on the KLM support the existence of a large FS. The height of the singularity in the 
electron momentum distribution decreases and becomes very small as the Kondo interaction 
goes towards the weak coupling region. These features are however very hard to observe. 
A more accurate study is certainly needed to confirm our conclusions. Finally, the nature 
of the heavy quasiparticles appear to be different from that proposed in the renormalized 
band structure studies [|r^. In the renormalized band picture, the heavy quasiparticles 
result from the small hybridization between the conduction electrons and the renormalized 
/—bands pinned at the Fermi level of the conduction electron sea. This requires the number 
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of /-electrons per site to be sligthly less than one. Our results suggest that these are instead 
loosely bound states made up of conduction electrons and /—spin fluctuations. This is 
consistent with recent results from field theory and exact diagonalization |T8| . 

We wish to thank S. Fujimoto for indicating ref. [6] to us. This work was supported by a 
grant from the Natural Sciences and Engineering research Council (NSERC) of Canada and 
the Fonds pour la formation de Chercheurs et d'Aide a la Recherche (FCAR) of the Quebec 
government. 
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FIGURES 

FIG. 1. The electron momentum distribution nc{k) for p — 0.75, 0.5 and 0.25 at Jk — 
and Jh — 0.5. 

FIG. 2. The magnetic structure factor Sf{k) for p = 0.75,0.5 and 0.25 at Jk = 10 
Jh = 0.5. 

FIG. 3. nc{k) for various values of at p = 0.5 and Jh = 0.5. 
FIG. 4. Sf{k) for various values of at p = 0.5 and Jh = 0.5. 
FIG. 5. Sf{k) for J^^ = 1.25, 1.5, 1.75, 2 and 2.5 at p = 0.75 and Jh = 0. 
FIG. 6. n^ik) for Jk = 1.5, 2 and 2.5 at p = 0.75 and Jh = 0. 
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